Superconducting density of states and vortex cores of 2H-NbS2 
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Scanning tunneling microscopy and spectroscopy (STM/S) measurements in the superconduct- 
ing dichalcogenide 2H-NbS2 show a peculiar superconducting density of states with two well defined 
features at 0.97 meV and 0.53 meV, located respectively above and below the value for the supercon- 
ducting gap expected from single band s-wave BCS model (A=1.76kBT" c =0.9 meV). Both features 
have a continuous temperature evolution and disappear at T c — 5.7 K. Moreover, we observe the 
hexagonal vortex lattice with radially symmetric vortices and a well developed localized state at the 
vortex cores. The sixfold star shape characteristic of the vortex lattice of the compound 2H-NbSe2 
is, together with the charge density wave order (CDW), absent in 2H-NbS2- 
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The study of the coexistence of superconductivity 
with competing physical phenomena such as magnetic or 
charge order has historically produced great interest on 
the scientific community. Anisotropics or modulations of 
the superconducting properties (in real and/or reciprocal 
space) often appear as a consequence of competing orders 
within the same system [TJ [2J. In the compound 2H- 
NbSe2, superconductivity appears within a CDW state 
(Tcdw=33K and T c =7.2 K) [3]. Low lying excitations 
measured deep in the superconducting state long time 
ago by specific heat [U [5] have been explained by recent 
experiments and theoretical calculations with a multi- 
band superconductivity and a peculiar anisotropy of the 
superconducting gap [6l [7] . Recent angular resolved pho- 
tocmission spectroscopy measurements demonstrate that 
the superconducting gap has, close to T c (at 5.7 K), 
largest values at k-space positions connected with CDW 
wavevectors [H [9]. Hess et al. [TU1 EEU H2] found that 
the local superconducting density of states (LDOS) at 
the center of the vortex core shows a high peak close to 
the Fermi level highlighting the lowest quasiparticle state 
bound within the vortex core well [13 . Around the vor- 
tex core, the LDOS is far from respecting in-plane sym- 
metry and intriguing vortex lattice images with patterns 
showing strong in-plane LDOS modulations are obtained 

2H-NbSe2 belongs to the transition-metal dichalco- 
genides (2H-MX 2 with M = Ta, Nb and X = Se, S), a 
family of systems which is unique to study the interplay 
between CDW order and superconductivity. The 2H- 
MX 2 compounds share a double layered structure made 
of two hexagonal X sheets with an intercalated M sheet 
(X-M-X), connected through very weak van der Waals 
bonds [14 . This produces highly anisotropic, quasi two 
dimensional electronic properties. The features of the 



Fermi Surface (FS) expected to be common in all sys- 
tems of the series are two concentric cylindrical FS sheets 
centered on both T and K points, derived from the 
transition-metal d bands [15j [16j HZl [18]. When going 
over the series from 2H-TaSe 2 and 2H-TaS 2 to 2H-NbSe 2 
and 2H-NbS 2 , the ratio of the intralayer lattice constant 
with the interlayer distance a/c increases, as well as T c , 
whereas Tcdw decreases strongly O [THl HO]- In 2H- 
NbS 2 , CDW order seems to be absent from macroscopic 
measurements and T c is around 6 K [T31 [5T] , very close to 
T c of 2H-NbSe 2 (7.2 K). The study and characterization 
of superconducting properties of 2H-NbS 2 arises as an ex- 
traordinary route to give new insight into the unsolved is- 
sues raised by the peculiar behavior of 2H-NbSe 2 . Among 
these issues we highlight the following: Is two band su- 
perconductivity common to both compounds? How do 
quasiparticles behave in the vortex core wells of both 
compounds? And finally (possibly the most intriguing 
question), what is, if any, the relationship between the 
star shaped vortex cores in 2H-NbSe 2 and CDW order? 
Here we discuss STM/S measurements of the supercon- 
ducting properties of 2H-NbS 2 which give relevant an- 
swers to these questions. 

We use a homebuild STM/S system installed in a di- 
lution refrigerator. The experimental setup, as described 
in previous work [551 1221 (HI H5j , is mechanically decou- 
pled of external vibrations, has a resolution in energy 
of 15 /zeV that has been previously tested by measur- 
ing clean spectra in low T c superconductors as such Al, 
and a mechanical pulling mechanism that allows to move 
the tip with respect to the sample holder. 2H-NbS 2 was 
synthesized as described elsewhere [5B]. Thousands of 
layered crystals with different shapes and sizes, typically 
of some tens of /Km 2 were obtained. The chemical com- 
position of these crystals was checked using a scanning 



electron microscope JEOL 5800 equipped with a micro- 
analyzer. From powder X Ray diffraction technique, the 
2H structure was found to be the dominant polytypic 
form of the batch and only a very small amount of 3R- 
NbS2 was detectable. Single crystal diffraction experi- 
ments performed at room temperature using a four cir- 
cles FR 590 Nonius CAD-4F Kappa-CCD diffractometer 
confirmed the presence of the 2H polytype and revealed 
that the platelets with a well defined hexagonal shape 
have less stacking faults along the c-axis than the others. 
For the present experiment, we have chosen three small 
single crystals with a nice hexagonal shape, which were 
individually mounted on a gold sample. It was possible 
to move easily the scanning window from sample to sam- 
ple crossing a clean surface of gold. The superconducting 
features discussed below are well reproducible in all three 
samples. 
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FIG. 1: a) Atomic scale topography (left panel), taken with 
STM in 2H-NbSe2 at O.IK. A superimposed modulation due 
to the presence of the CDW is observed. The corresponding 
Fourier transform image (right panel) shows the wave vectors 
of the atomic and CDW modulations in the reciprocal space 
(marked by arrows in the inset), b) Atomic scale topography 
at 0.1 K (left panel) and its Fourier transform image (right 
panel) at atomic scales in 2H — NbS2- No CDW order is 
observed. 



STM images in 2H-NbS2 show atomically flat regions 
with lateral sizes that easily go up to several hundreds 
of nanometers. In fig. [I] we present a STM image with 
atomic resolution obtained at 0.1 K in 2H-NbS2 (fig. 
[lb), and compare it to a STM image in 2H-NbSe2 mea- 
sured under the same conditions (fig. [Ik). In 2H-NbSe 2 , 
we observe the hexagonal surface atomic lattice of Sc 
atoms and the CDW order as previously reported in 
many other STM experiments and largely discussed in 
literature [27l[28]. In 2H-NbS 2 , the STM image (fig. [lb) 
also shows the hexagonal symmetry of the surface atomic 
lattice of the chalcogenide (here S), characteristic of this 
family of materials, with a slightly smaller lattice con- 



stant, as expected from X-Ray measurements [2"Tl |2"61, |2"9"] . 
However, both the STM image and its associated Fourier 
transform in 2H-NbS2 (left and right panels in fig. [lb, re- 
spectively) do not show any modulations besides of the 
ones coming from atomic periodicity. This result gives 
strong evidence that there is no coexistence of supercon- 
ducting and CDW orders in 2H-NbS2 at least down to 
0.1 K. 



o 
7 




-2-1012 
Bias voltage (mV) 



FIG. 2: Characteristic tunneling conductance curve obtained 
in 2H-NbS2 and measured at O.IK. The derivative of the as- 
sociated LDOS (shown in the inset) presents two peaks, at 
0.97meV and 0.53meV (marked by blue and green arrows, 
respectively), which correspond to two peaks in a gap distri- 
bution. 



In fig. [2] we show a typical tunneling conductance vs. 
bias voltage curve obtained in an atomically flat region at 
0.1 K. The curve is fully reproducible over the whole sur- 
face, showing tiny differences (below 10 %) as a function 
of the position at atomic scale, in a similar way as found 
in 2H-NbSe2 previously [21]. As is well known, the lo- 
cal tunneling conductance simply reflects the convolution 
between the LDOS of the sample and the derivative of 
the Fermi function. At these temperatures, the tunnel- 
ing conductance curve pretty directly follows the LDOS, 
because the derivative of the Fermi function is very nar- 
rowly peaked, in particular as compared to the critical 
temperature (here 0.1K s=s 0.02T C ). The LDOS is very far 
from the divergent quasiparticle peaks expected within 
single band s-wave BCS model, and shows instead a slight 
increase from relatively low energies, of about 0.3 meV, 
which becomes much more marked around 0.6 meV, and 
a rounded quasiparticle peak somewhat above 1 meV. 
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FIG. 3: a) Temperature scan of tunneling conductance spec- 
tra (left) and superconducting LDOS (right). The LDOS 
is obtained by deconvolution from tunneling conductance 
curves. The data are, from bottom to top, taken at 0.1K, 
0.9K, 1.5K, 2K, 2.5K, 3K, 3.5K, 4K, 4.25K, 4.5K, 4.75K, 5K, 
5.2K, 5.4K and 5.6K. Red lines in left panel are the fits ob- 
tained from convoluting the LDOS shown in the right panel 
with temperature, b) Temperature dependence of the two 
peaks in d(LDOS)/dE (green and blue points, corresponding 
to the features marked by green and blue arrows in fig. [2] 
dashed lines are guides to the eye). The solid line is the BCS 
expression taking A = 1.76k B T c = 0.87 meV with T c = 5.7 
K. 



Such a LDOS evidences that more than a single gap is 
open over the Fermi surface, i.e. there is distribution of 
values of the superconducting gap. It is easily shown that 
d(LDOS)/dE has peaks at energies which correspond to 
the superconducting gaps which appear more often in 
the distribution. As shown in the inset of fig. [2] two 
well defined peaks d(LDOS)/dE are visible at energies 
±(0.97 ± 0.03) meV and ±(0.53 ± 0.03) meV. The tem- 
perature dependence of the tunneling conductance and 
of the LDOS (obtained as in Ref. [23]) is shown in fig. 
[3^. The LDOS has a smooth temperature dependence. 
The features appearing at low temperatures are simply 
scaled to lower energies as T c is approached. The tem- 
perature dependence of the two peaks in d(LDOS)/dE 
is shown in fig. [3b. Both of them disappear at T c =5.7 



K. Moreover, the value of the gap expected within single 
gap s-wave BCS model (1.76 kBT c =0.87 meV; the solid 
line in fig. [2b gives its temperature dependence), is in be- 
tween the values of the two peaks on d(LDOS)/dE at all 
temperatures (fig. [2b). This gives strong evidence that 
two different superconducting gaps (±(0.97 ± 0.03) meV 
and ±(0.53 ± 0.03) meV) open at different parts of the 
Fermi surface of 2H-NbS2- 

The situation resembles the results found previously 
in MgB 2 [3D1 HU |33 |33 |31 |35], and in 2H-NbSe 2 
[3 HI 13 123 H3 [35], which evidenced the appear- 
ance of two band superconductivity in these materials. 
In 2H-NbS2 the gap values appear to be much closer to 
each other than in the case of MgB 2 , where the two gaps 
differ by more than a factor of three. Moreover, inter- 
band scattering effects should be present to allow for a 
smooth temperature dependence of both gaps [37] . How- 
ever, in 2H-NbS2 the two gap nature seems better de- 
veloped than in 2H-NbSe2, where similar features in the 
density of states are less clearly resolved [7J [TTJ HO] ■ 
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FIG. 4: a) STS image of the vortex lattice in 2H-NbS 2 taken 
at 0.1 K and 0.15 T is shown in a (360 nm x 360 nm). In b, we 
show, for comparison, a STS image of the same area obtained 
in 2H-NbSe2 under similar conditions. In c we present tun- 
neling conductance curves in 2H-NbS2 along a line of about 
65 nm that goes from the middle point between vortices into 
the vortex core (using a similar color scheme as in a and b), 
following the line shown in a. 

Under magnetic field we observe an ordered hexagonal 



vortex lattice in atomically flat regions for different val- 
ues of the applied field. A STS image of the vortex lattice 
at 0.1 K and 0.15 T built from the conductance curves at 
zero bias voltage is shown in the fig. HK, and compared 
to results at the same field and temperature in 2H-NbSe 2 
(fig. [4]d) rniJUniH]. The vortex core shapes, as observed 
with STS, are dramatically different in both compounds. 
First, the star shape of 2H-NbSe2 is fully absent in 2H- 
NbS 2 . When making STS images at different energies, 
vortex cores in 2H-NbS2 appear to be round, fully re- 
specting angular symmetry, at all energies, whereas 2H- 
NbSe2 has a rich phenomenology, which includes the for- 
mation of complex angular patterns, discussed in detail 
in previous works fTUJ [TTJ [12] . On the other hand, the en- 
ergy dependence of the LDOS at the center of the vortex 
core in 2H-NbS2 is very similar to the one found in 2H- 
NbSe2- There is a high peak at the Fermi level, due to 
the lowest level of bound localized states inside the vortex 
core wells [13], which splits when going out to the vor- 
tex core, exactly as in 2H-NbSe2. The difference between 
both compounds lies only in the angular dependence of 
these features in the LDOS. Whereas in 2H-NbSe2, the 
way the localized state peak splits is complex and shows 
a very strong real space angular dependence, in 2H-NbS2, 
the LDOS respects in-plane symmetry. This must be cor- 
related to the absence of CDW in 2H-NbS 2 (Fig.l), so we 
can state that the peculiar star shape of the Abrikosov 
vortex lattice in 2H-NbSe2 and its voltage dependence is 
due to the presence of CDW order in this compound. 

Note that two band superconductivity appears to be 
common to 2H-NbSe2 and to 2H-NbS2, and it could be 
a rather general feature of this family of materials which 
deserves further study. For example, two different gaps 
may open at the two groups of transition metal d-electron 
derived FS sheets[18]. However, 2H-NbSe2 has an addi- 
tional gap modulation caused by CDW that shapes the 
angular dependent properties of the different supercon- 
ducting gaps opened over one of the bands, or eventually 
over the whole Fermi surface. Indeed, different models 
directly associate the star shape in the vortex core local- 
ized state features in 2H-NbSe2 with a strong in-plane 
modulation of the superconducting gap [351 IM1 HO] • The 
superconducting properties of 2H-TaS2 and of 2H-TaSe2, 
where CDW order appears at a much higher tempera- 
ture (80 K and 120 K respectively), could be even more 
intriguing and may also include a richly shaped vortex 
lattice. 

On the other hand, naturally, the presence of features 
in the LDOS due to localized bound states strongly de- 
pends on the mean free path, as scattering flattens out 
the peak in the LDOS at the center of the vortex core 
[411 142] , So the observation of the peak in the LDOS 
in 2H-NbS2 at the center of the cores shows that our 
samples are within or close to the clean limit. Also, 2H- 
NbS2 appears as the system that best approaches the 
initial proposal of Caroli, de Gennes and Matricon re- 



garding quasiparticle bound states inside cores [T3]. In- 
deed, whereas 2H-NbSe2 is strongly influenced by CDW, 
vortex core features in other systems, as e.g. high T c 
cuprate and related superconductors, boron doped di- 
amond, or YM2B2C have a rather rich and intriguing 
behavior [131 3U US] , with, however, no clear bias sym- 
metric peak appearing at the Fermi level LDOS at the 
center of the cores. 

In summary, we have found that the LDOS of 2H- 
NbS2 has two pronounced features related to intrinsic two 
gap superconductivity, in close analogy to MgB 2 . Vortex 
cores show a bound state close to the Fermi level and 
do not present any in-plane anisotropy. By comparing 
2H-NbS2 with previous work on 2H-NbSe2, we find that 
charge order is the main responsible of the strong in- 
plane six-fold anisotropy that occurs in 2H-NbSe 2 . In 
particular, CDW is responsible for the six fold structure 
of the vortex lattice found by Hess et al. [TT] in 2H-NbSe2 . 
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